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Abstract We consider n-dimensional discrete motions such that any two neighbouring positions correspond 
in a pure rotation ("rotating motions"). In the Study quadric model of Euclidean displacements these motions 
correspond to quadrilateral nets with edges contained in the Study quadric ("rotation nets"). The main focus 
of our investigation lies on the relation between rotation nets and discrete principal contact element nets. 
We show that every principal contact element net occurs in infinitely many ways as trajectory of a discrete 
rotating motion (a discrete gliding motion on the underlying surface). Moreover, we construct discrete 
rotating motions with two non-parallel principal contact element net trajectories. Rotation nets with this 
property can be consistently extended to higher dimensions. 
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1 Introduction 

Discrete differential geometry is an active field of ge- 
ometrical research. Its aim is the development of dis- 
crete notions for well-known concepts from differential 
geometry. The resulting theories are often more ele- 
mentary and concrete when compared to their smooth 
counterparts. While classic differential geometry is 
largely based on analysis, elementary geometric inci- 
dence or closure theorems are at the core of discrete 
differential geometry. It is a discipline that naturally 
lends itself to applications that require numeric sim- 
ulation, visualization, or the building of real world 
objects. An excellent introduction to the current state 
of research is the monograph 0. 

In this article we relate recent progress in the the- 
ory of discrete curvature line parametrizations to spa- 
tial kinematics. We study discrete nets of proper 
(orientation-preserving) Euclidean displacements such 



that neighbouring positions correspond in a relative ro- 
tation ("discrete rotating motions"). Smooth motions 
with this property naturally arise as gliding motions 
along principally parametrized surfaces. Their discrete 
counterparts, rotating motions with discrete curvature 
line trajectories, are the main topic of this article. 

In Section [2] we recall the notion of principal con- 
tact element nets — families of contact elements (point 
plus oriented tangent plane), indexed by Z", such that 
neighbouring contact elements have a common tangent 
sphere. Principal contact element nets have been intro- 
duced in as a comprehensive concept that captures 
different notions of discrete principal parametrizations 
(circular nets and conical nets, see [2, Section 3.1] and 
flU). Indeed, the points of a principal contact element 
net form a circular net (the elementary quadrilaterals 
are circular; Figure [TJ left), while its planes form a 
conical net (four planes meeting in a vertex are tangent 
to a cone of revolution; Figure [TJ center). 
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A concise formulation of all calculations and for- 
mulas in this article is possible within the dual quater- 
nion calculus of spatial kinematics (Section 2.2 1. In 
Section 1231 we recall fundamental results on "rotation 



quadrilaterals" [ 10 ], the elementary building blocks of 
rotation nets. 

The major contributions of this article are presented 
in Section[3] Just as a smooth surface gives rise to many 
gliding motions, parametrized by principal lines and 
the rotation angle about surface normals, a principal 
contact element net occurs in many ways as trajectory 
of a discrete rotating motion. The main result is a 
proof of the multidimensional consistency of discrete 
rotating motions with two independent principal con- 
tact element trajectories. The two trajectory surface 
are related by discrete version of the classic Backlund 
transform for pseudospherical surfaces. 

2 Preliminaries 

2.1 Curvature line discretizations 

A parametrization of a smooth surface in M 3 is 
called a curvature line parametrization or principal 
parametrization if infinitesimally neighbouring sur- 
face normals along both families of parameter lines 
intersect. Generically, this parametrization is unique 
(up to re-parametrization of the individual parameter 
lines). The property of concurrent neighbouring nor- 
mals is preserved in its usual discretizations. The most 
prominent example of discrete curvature lines are circu- 
lar nets — quadrilateral nets such that any elementary 
quadrilateral has a circumcircle, see for example J2J 
Section 3.1]. An alternative discretization are coni- 
cal nets — quadrilateral nets such that the planes meet- 
ing in a vertex are tangent to a cone of revolution (H. 
Neighboring circle axis and neighbouring cone axis 
in circular and conical nets intersect and can serve as 
discrete surface normals. In Q principal contact ele- 
ment nets were introduced as a generalization of both, 
circular and conical nets. 

Definition 1. An oriented contact element is a pair 
(jj, v) consisting of a point p and an oriented plane v 
incident with p. The oriented line N orthogonal to v 



and incident with p is called the axis or normal of the 
contact element, p is its vertex and v its tangent plane. 

Contact element nets are quadrilateral nets of ori- 
ented contact elements, that is, they are maps from Z" 
to the space of oriented contact elements. The image 
of i 6 Z" is denoted by (p-„ Vi). In Definition [2] below 
we adopt the notation of (2] where T, indicates a shift 
of indices in the i-th coordinate direction. For example 
= P220, T2P120 = P130, etc. 

Definition 2. A contact element net is a map i h-> 
(pi,v{) from W to the space of oriented contact ele- 
ments. A principal contact element net (Figure [TJ is 
a contact element net such that any two neighbouring 
contact elements (p-„ Vj), ^i(pi, Vj) have a common ori- 
ented tangent sphere. 

In other words, neighbouring normals N-, and TjiVj 
in a principal contact element net intersect in a point 
z[ = N; n tiNi which is at the same oriented distance 
from pi and Z\p\. The points p\ in a principal contact 
element net constitute the vertices of a circular net, 
the oriented planes Vj are the face planes of a conical 
net. The contact element axes form a discrete line 
congruence (21 Section 2.2] such that any elementary 
quadrilateral lies on a hyperboloid of revolution (Fig- 
ure^ right). Discrete line congruences of this type 
have interesting properties but have not yet been dis- 
cussed in literature. 

In this article we consider a kinematic generation of 
principal contact element nets in R 3 . A concise ana- 
lytic description can be obtained by means of the dual 
quaternion calculus of spatial kinematics which shall 
be introduced now. 

2.2 Kinematics and dual quaternions 

Quaternions and dual quaternions are important tools 
in theoretical and applied kinematics, see for example 
the description in Q, J6j Section 4.5] or iTTTl Chap- 
ter 9]. In our study they turn out to be a versatile tool 
as well. Equations in dual quaternion form are concise, 
of low degree, accessible to geometric interpretations 
and free of the need of case distinctions. We assume 
that the reader is familiar with quaternion algebra (see 
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Figure 1: Elementary quadrilateral in a principal contact element net: Points on a circle (left), planes tangent to 
a cone of revolution (middle), normal lines on a hyperboloid of revolution (right) 



for example or ||6j Chapter 4]) and describe only 
the extension to dual quaternions and its application to 
spatial kinematics. 

A dual quaternion is an object of the form a = a + ed 
where primal part a and scalar part a are ordinary 
quaternions and e is the dual unit satisfying e 2 = 0. The 
addition of dual quaternions is performed component- 
wise for primal and dual parts. The dual quaternion 
multiplication * extends the multiplication of ordi- 
nary quaternions. It is associative, distributive and 
the quaternion units 1, i, j, k commute with e. These 
properties define the dual quaternion multiplication 
uniquely. 

A dual quaternion a = a + ed can be identified with a 
vector a = (a, a) = (ao, ■■■ , a-]) in IR 8 . Vector part and 
dual part of a are 

Q3a = (0,ai,a2,aj,0,a5,a(,,aj) = QJa + gQJa, 

ea = (a , 0,0,0,04,0,0,0) = Sd+eGd. ' ' ' 

A normalized dual quaternion a = a + ed satisfies two 
conditions: 

\\d\\ 2 = a} j + a\ + a\ + a\ = \ and (a,a) = (2) 
where 

3 

(a,b) := £(a,-6i + 4 + a i4 4&i). (3) 

i=0 



The second equation in ([2]) is the well-known Study 
condition. 

Plucker coordinate vectors P = (po, . . . ,p$) for 
straight lines (see for example |7, Chapter 2]) can be 
embedded in the space of dual quaternions as 

P^->(0,po,pi,P2,0,P3,P4,P5)- (4) 

In this case the Study condition Q reduces to the 
Plucker condition such that normalized Plucker coor- 
dinates (characterized by p\ + p\ + p\ = 1) become 
normalized dual quaternions of vanishing scalar part. 
In this article we do not distinguish between a straight 
line and its Plucker coordinates, embedded in the space 
of dual quaternions. 

Normalized dual quaternions with non-vanishing pri- 
mal part constitute a two-fold covering of the group 
of proper Euclidean displacements. The action of a 
normalized dual quaternion a on a point with coordi- 
nate vector (jci,X2,X3) and a line P (both in the moving 
space) is given by 

1 + ex' = a e *(\ + £x)*a~ l , 

n -i ^ 
P = a e *P*a e , 

where x = (0,X\,X2,X3) is a vector valued ordinary 
quaternion, a e :=a — ea, and the prime indicates coor- 
dinates in the fixed space. The inverse displacement is 
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described by the inverse dual quaternion a denned 
by the equation a-kaT x = 1. It is uniquely defined for 
all dual quaternions with non- vanishing primal part. 

Since a and —a describe the same displacement, 
it is natural to identify proportional dual-quaternions, 
thus arriving at Study's kinematic mapping which asso- 
ciates proper Euclidean displacements with points of 
the Study quadric 

y-. (x,x)=0. (6) 

The Study quadric is a hyperquadric in the seven- 
dimensional projective space P 1 over M 8 . Only points 
of an exceptional three-space E with equation x = 
(0,0,0,0) do not occur as images of proper Euclidean 
displacements. We always identify Euclidean displace- 
ments with homogeneous coordinate vectors that de- 
scribe points on the Study quadric (minus E). 

2.3 Rotation quadrilaterals 

In Section [3] we will consider special quadrilateral nets 
in the Study quadric. The geometry of an elementary 
quadrilateral 

[a^tia^Tja^TiTja^ = [a ,ai,a 2 ,a 3 ) (7) 

of a net of this type is discussed in this section. It 
follows the presentation of [ 10]. 

Definition 3. A quadrilateral [ao,a\ ,02,03] on the 
Study quadric is called a rotation quadrilateral if its 
vertices and edges are contained in y . 

The name "rotation quadrilateral" is justified by the 
observation that the edge through a; and is con- 
tained in y if and only if the relative displacement 
*"(',;+ 1 := *a7 is a pure rotation (indices modulo 
four; see JT2l Satz 19] or Q). (We will frequently use 
results and formulas of lPT2l but, occasionally, adapt 
them to match our convention of quaternion multiplica- 
tion which is slightly different from that used in [ 12].) 
The algebraic characterization of relative rotations is 

X5 (ru+i) = X5 (fli+l * aj 1 ) = (8) 

where 715 (x) is the projection onto the fifth coordinate 
(the dual scalar part) of a dual quaternion [12, Satz 13]. 



We denote the relative revolute axis of r, j + \ in the mov- 
ing space by R^+i. 

The main result of iTTOll characterizes contact ele- 
ments whose homologous images form an elementary 
quadrilateral in a principal contact element net. 

Proposition 4 ( H10I1 ). The only contact elements in 
the moving space whose homologous images with re- 
spect to a generic rotation quadrilateral form a non- 
degenerate elementary quadrilateral of a principal con- 
tact element net are those whose axes are transversal to 
the four relative revolute axes Rq\, Rn, R23> and R30. 

Further elementary quadrilaterals of principal con- 
tact element nets have the normal Rij+i. Since their at 
and a, + i images are identical, these quadrilaterals are 
degenerate. For a generic rotation quadrilateral there 
exist two (possibly complex or coinciding) transversals 
M and N of the four relative revolute axis. An impor- 
tant property of rotation quadrilaterals is stated in the 
following completion theorem: 

Theorem 5. Consider two skew lines M, N in the mov- 
ing space and two displacements ao, a 2 of a rotation 
quadrilateral. Generically, there exist two positions 
a\, c?3 (possibly complex) such that [ao,a\,a2,as] is a 
rotation quadrilateral with relative revolute axes Ra+\ 
that intersect M and N. 

Proof. We consider a (yet undetermined) position x 
such that the relative displacements x^a^ 1 , x-ka^ 1 are 
rotations whose axes Xo, X2 intersects M and Af. It turns 
out that this problem amounts to solving a quadratic 
equation from which the Theorem's claim follows. 

The Plucker line coordinate vector of the relative 
rotation axis X- in the fixed space is 

Xi = *8(xi) e (9) 

(adapted from H2J Satz 13]). According to Q, the 
Plucker coordinate vector of the relative revolute axis 
Xi in the moving space is 

Xi = (a i )- 1 *Xl*(a i )e. (10) 

Hence, the sought position x has to satisfy the six linear 
equations 

{{a i ) e -k^O{x-kar l ) e -k(a i )-\T) = 0, 
7is(x-kaJ l ) = 
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with i G {0, 2} and T G {M, Af}, and the quadratic equa- 
tion (x,x) = 0. The solutions are the intersection points 
of a straight line with the Study quadric ([6]). □ 

We will actually need the result of Theorem[5]in the 
following form: 

Corollary 6. It is possible to complete a rotation 
quadrilateral whose relative revolute axes intersect 
two skew lines M, N in the moving space from three 
admissible input positions ao, a\, a%. The input data 
is admissible if the two relative displacements a\-ka~^} 
and a2*a^[ l are rotations whose axes intersect M and 
N. In this case, the missing position a^ is unique and 
real (if the input positions are real). 

3 Discrete rotating motions 

Now we are ready to introduce the central concept of 
this article: 

Definition 7. A rotation net (or a discrete rotating mo- 
tion) is a quadrilateral net whose vertices and edges are 
contained in the Study quadric. 

The elementary quadrilaterals of rotation nets are 
rotation quadrilaterals. Rotation nets of dimension 
two discretize two-parameter motions x{t\,t2) with pa- 
rameter lines tj = const, whose instantaneous screws 
have vanishing pitch, that is, they are actually rota- 
tions. The geometric interpretation in terms of the 
Study quadric is that not only the point x(t\ , tj) but also 
the tangents to the parameter lines are contained in the 
Study quadric. We call parametrized motions of this 
type rotating as well. The extension of this concept to 
more-dimensional motions is obvious. 

Important examples of two-dimensional rotating mo- 
tions are obtained from a principal parametrization 
f(?i,?2) of a surface <J> C M 3 . The Darboux frame asso- 
ciated with this parametrization is the trihedron with 
base f and legs ui, 112, and n where 



df/dtj 

W/dti 



n = ui x 112. 



(12) 



The Darboux frame is orthonormal. As t\ and ti vary 
in their respective domains, it moves along the surface 



<£. The thus defined motion x{t\,t2) is rotating and 
we may call it the Darboux motion of the principal 
parametrization. If we add the rotation angle w about 
the surface normal n(?i ,t-i) as a third parameter the mo- 
tion x(t\,t2,w) is still rotating (see J7J Section 7.1.5], in 
particular Remark 7.1.16). It is called a gliding motion 
along <!>. Similar results hold true for more-parameter 
motions. These observations motivate our study of the 
relation between discrete rotating motions and discrete 
principal curvature parametrizations. 

3.1 Principal contact element nets and rotating 
motions 

We are going to discuss possibilities for generating 
a principal contact element net (pi,Vi) as trajectory 
of a discrete rotating motion. Any two neighbouring 
contact elements (pi, Vj) and T/(/?|, Vj) have a plane of 
symmetry j6; and, starting from one contact element, 
the complete net can be generated by successive re- 
flections in the planes /$; . The sequence of reflections 
between two contact elements is not unique but never 
leads to contradictions. 

Sometimes (as for example in 0) it is useful to as- 
sign an orthonormal frame (Xi,Yi,Zi) to every element 
of a principal contact element net such that the ori- 
gin coincides with p\, Z\ is the normal line and the X,- 
and Fj-axes correspond in the reflections at the planes 
j3;'. The thus obtained discrete line congruences X\ and 
Y\ can be regarded as discretizations of one family of 
principal curvature lines. It is obvious that by a sim- 
ple change of orientation of certain X r and Fj-axes the 
reflection at j6; can be replaced by a rotation about 
an axis perpendicular to N\ and ZiNi and incident with 
z! = Ni D TjNi G /3/. The thus obtained rotation net can 
be considered as discretization of a Darboux motion. If 
we perturb every frame by a rotation through a certain 
angle about its Zj-axis, neighbouring frames still corre- 
spond in a rotation about an axis through z,[ and in jS/. 
This rotation net can be considered as a ^-dimensional 
discrete gliding motion along a principal parametriza- 
tion. 

These considerations show that every principal con- 
tact element net occurs in multiple ways as trajectory 
of a rotation net — just as any principal parametriza- 
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tion gives rise to infinitely many gliding motions. The 
rotation angle about the Zj-axis gives one degree of 
freedom per vertex. 

3.2 Pairs of principal contact element nets 

By Proposition [4j rotation nets with principal contact 
element nets as trajectories are characterized by the 
fact that all relative rotation axis in the moving space 
intersect a fixed line N whose images constitute the 
set of contact element normals. Clearly, every point 
p £ N and the plane v through p and orthogonal to 
N define a contact element (p,v) whose trajectory is 
a principal contact element as well. In other words, 
principal contact element nets as trajectories of discrete 
rotating motions come in one -parameter families of 
parallel nets. 

Calling two principal contact element nets indepen- 
dent if there exists a pair of contact elements with the 
same index i £ Z" but different normals, we aim at a 
kinematic generation of independent principal contact 
element nets as trajectories of rotation nets. In view of 
Proposition [4] we can hope at most for two independent 
trajectories. The characteristic property is that all rel- 
ative rotation axes in the moving space intersect two 
fixed lines M and N. We will show that such rotation 
nets exist for arbitrary dimension of the underlying 
motion. This result is rather surprising since a naive 
counting of free parameters suggest only existence of 
2-dimensional motions with this property. 

Theorem 8. A discrete n-dimensional rotating motion 
with two independent principal contact element nets 
as trajectories is uniquely defined by two skew axes M, 
N of the contact elements in the moving space and the 
values along the coordinate axes in W. 

For n = 2 this result follows directly from Corol- 
lary [6] Starting with aoo, oio, and aoi, the position an 
is uniquely defined. From aio, 020, and a\\ we can 
construct az\, and so on. Two independent principal 
contact element net trajectories of a discrete rotating 
motion of dimension n = 2 are depicted in Figure [2] 
The figures formed by corresponding points and face 
normals are congruent. Neighbouring figures corre- 
spond in a pure rotation (Figure [3]). 



In case of n> 3 it is not immediately clear that this 
inductive construction works. We describe the situa- 
tion only for n = 3; the problems in higher dimensions 
are similar. Consider the elementary cube a^ with 
i, j, k£ {0, 1} (Figure]?]). The input data consists of 
oooo, « ioo, ao\o, and c?ooi- According to Corollary [6] it 
defines ano, oioi, and aon- Now there are three possi- 
bilities to construct the missing vertex am, from the 
three positions aioo, ano, aioi, from the three positions 
ooio, ai io, aon, or from the three positions aooi, oioi, 
flon- We have to show that all thus constructed posi- 
tions are actually identical. In the terminology of 
this is called 3D -consistency of rotation nets with two 
independent principal curvature trajectories. In Theo- 
rem |9| below, we will show that these nets are actually 
nD-consistent. This is a fundamental property in the 
discretization of differential geometric concepts and 
immediately implies Theorem [8] 

Theorem 9. Generic discrete rotation nets with 
two independent curvature line trajectories are nD- 
consistent. 

Proof. We have to show that a generic ^-dimensional 
cube can be constructed from one vertex ar Q ) and 
^-adjacent vertices T/ajo o)> i € {1, . • • ,«}■ Using the 
notation of the last paragraph before the Theorem we 
proof 3D-consistency at first. We already know that 
there exist points a\, ai and a 3 such that 

[aioo,ano,aioi,fli], 
[0010,0110,0011,02], (13) 
[aooi, aiobooii, 03] 

are rotation quadrilaterals with all required properties. 
We have to show that these points coincide whereupon 
we may set sm := a\ = c?2 = 03. 

If we require only one principal contact element 
net trajectory (for example with normal M and plane 
/J. _L M), only six linear equations remain. They define 
a straight line K C P 1 which, in the sense of algebraic 
geometry, has at least two intersection points k\, &2 
with the Study quadric. Since principal contact element 
nets are known to be «D-consistent [3], the k{- and 
&2-images of the contact element (p,n) coincide. We 
denote this contact element by (pk,I~Lk) and conclude 
that the relative displacement between k\ and &2 is a 
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pure rotation. This implies that the line K is actually 
contained in the Study quadric 5?. 

Denote by L C ^ the straight line of positions ob- 
tained in the same way but with M replaced by N and by 
(Pl,Vl) the corresponding contact element. The proof 
of 3D-consistency will be finished if we can show that 
K and L have a point in common. Assume conversely 
that K and L are skew. Then they span a three- space 
whose intersection with 5? is a ruled quadric. We con- 
clude that every position k G K can be rotated to a 
unique position / G L. In other words, the contact ele- 
ment (pk , Vk) can be rotated in infinitely many ways 
into the contact element (pi, Vl). By elementary geo- 
metric reasoning this is only possible if both contact 
elements have a common tangent sphere. This contra- 
dicts the skewness of M and N. Hence, k and / intersect 
in a unique position am which implies the Theorem's 
statement for n = 3. 

Now we consider 4D-consistency (see Figure [5] 
which displays the projection of a 4D-cube). According 
to Corollary [6] the input data 

«oooo, tfiooo, ^oioO; aooiO) floooi (14) 
defines the positions 

flnoo, «ioio, aiooi, tfono, aoioij a oon- (15) 

By the previous discussion, we can construct the posi- 
tion 

aniO) anoi) tfioii, «oin (16) 

without a contradiction. Now we have four ways 
to construct the missing position aim by complet- 
ing a 3D-cube. In Figure [5] two of these cubes are 
highlighted. They share the common quadrilateral 
[aoi 11,^001 1,^1011, ami]- Thus, by Corollary [6) com- 
pleting the 3D-cubes leads to the same position for 
flu 11. This situation does not change if we consider 
other pairs of 3D-cubes. 

The same argument yield TiD-consistency: Induc- 
tively it can be shown that the input data defines all 
positions uniquely with exception of one position a\. 
This position can be constructed by completing n cubes 
of dimension n — 1. But any two of these cubes share 
a face of dimension n — 2 > 2 which already uniquely 
determines a\. □ 



The construction of ^-dimensional rotation nets as 
in Theorem [8] from the given input data is inductive. 
At every step it requires solving the equation system 
( [TT] ) augmented with the Study condition Q. The solu- 
tion is generically unique and can be computed linearly. 
When prescribing the positions on the coordinate axes 
only a subset of ( [TT] ) needs to be fulfilled. 

4 Conclusion and future research 

In an attempt to relate recent progress in the field of 
discrete differential geometry to kinematics this arti- 
cle introduces discrete rotating motions and studies 
possibilities to obtain curvature line discretizations as 
trajectories. 

We are already in a position to present some of the 
implications of our study. Assume that (xi,V\) and 
are independent principal contact elements ob- 
tained as trajectory surfaces of a discrete rotating mo- 
tion. The figures formed by corresponding contact 
elements are congruent (Figure [3} and clearly we can 
find two families of parallel trajectory surfaces with the 
same properties. Distinguished examples are obtained 
by choosing the contact points x and y on the common 
perpendicular of the contact element normals M and ,/V 
in the moving space. In this case x\ G pL\ and y\ G Vi and 
the trajectory surfaces satisfy all geometric properties 
of the classic Backlund transform for pseudospherical 
surfaces (surfaces of constant Gaussian curvature): 

• Corresponding points are at constant distance (inde- 
pendent of i G Z"), 

• corresponding tangent planes intersect in a constant 
angle (independent of i G Z M ), and 

• the connecting line of corresponding points lies in 
both tangent planes. 

This observation leads to the conjecture that the tra- 
jectory surfaces are discrete surfaces of constant Gaus- 
sian curvature in the sense of UJ. Indeed, we have a 
proof for this which shall be published elsewhere. By 
a result of UJ, this implies that the parallel trajectory 
surfaces are linear Weingarten surfaces. An analytic 
description of a Backlund transform but for a differ- 
ent discrete curvature line parametrization ("discrete 
O-surfaces") can be found in [9]. Our main results, the 
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existence of discrete rotating motions with two indepen- 
dent principal contact element net trajectories and their 
multidimensional consistency, are meant to provide 
the basis for the description of Backlund transforms of 
pseudospherical principal contact element nets — in a 
forthcoming publication. 
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